Resonance wave functions located at the Stark saddle point 
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We calculate quantum mechanically exact wave functions of resonances in spectra of the hydro- 
gen atom in crossed external fields and prove the existence of long-lived decaying quantum states 
localized at the Stark saddle point. A spectrum of ground and excited states reproducing the nodal 
patterns expected from simple quadratic and cubic expansions of the potential in the vicinity of the 
saddle point can be identified. The results demonstrate the presence of resonances in the vicinity of 
the saddle predicted by simple approximations. 

PACS numbers: 32.60.-K, 32.80.Fb, 82.20.Db 



I. INTRODUCTION 

The hydrogen atom in crossed static electric and mag- 
netic fields is an important example for a quantum sys- 
tem accessible both in experiments and numerical calcu- 
lations. Adding external fields to the Coulomb poten- 
tial of the hydrogen atom opens the possibility for wave 
functions to be localized far away from the nucleus. The 
existence of such localized states in quantum systems has 
attracted the attention of both theoretical [lMl] and ex- 
perimental 0,0] investigations over a long period of time. 
First investigations considered a gauge dependent variant 
of the system, in which the paramagnetic term was omit- 
ted [H 0, Q due to its smallness as compared with the 
diamagnetic term at large field strengths. The resulting 
sum potential possesses a minimum away from the nu- 
cleus in which localized states can exist. However, if one 
takes into account all terms of the magnetic field, one 
can immediately see that the outer potential minimum 
does not exist in the real physical system, rather, one 
finds classical electron motion close to the Stark saddle 
point which is bound in two directions and unbound in 
the third direction. In fact, localized resonance quan- 
tum states in the vicinity of the saddle point have been 
predicted by Clark et al. [|[ as quantized version of this 
quasi-bound motion. 

The question whether or not resonances located at the 
Stark saddle point do exist has also a crucial meaning 
for the ionization mechanism. Classical electron paths 
describing an ionization must pass the vicinity of the 
Stark saddle. This process was investigated, e.g., using 
the transition state theory 04111 • The transition state 
theory is applicable to many dynamical systems which 
evolve from an initial to a final state and has its origin in 
the calculation of chemical reaction rates. The concept is 
based on classical trajectories which describe a reaction 
by passing a boundary in phase space, viz. the transition 
state, which must be crossed by all trajectories connect- 
ing the initial ( "reactants" ) and final ("products") side. 
For the hydrogen atom in crossed fields the theory postu- 
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lates classical orbits confined in the vicinity of the Stark 
saddle point Q. 

The work by Clark et al. [H can be regarded as a first 
step in discussing the ionization mechanism of the hydro- 
gen atom in crossed electric and magnetic fields in this 
context. The quasi-bound states confined to the vicin- 
ity of the saddle found by Clark et al. coincide with an 
approximation to the Hamiltonian in the framework of 
the transition state theory, however, only an algorith- 
mic procedure based on a normal-form representation of 
a power-series expansion of the Hamiltonian allowed for 
identifying the transition state @, OH [H| ■ Even though 
the transition state for the hydrogen atom in crossed 
fields has been found, the question of whether or not 
the classical trajectories in its vicinity leave a signature 
in the exact quantum spectrum remained unanswered for 
a long time. Finally, clear evidence for signatures in the 
energies, i.e., the eigenvalues of the quantum resonances, 
was found jl2| . 

It is the purpose of this paper to demonstrate that 
there is, indeed, a relationship between the exact quan- 
tum resonances and the quantized energy levels of the 
electron motion near the transition state due to a defi- 
nite spatial restriction of the resonances' probability den- 
sity to a small region around the saddle. We calculate 
the position space representation of the corresponding 
wave functions and are able to show that they are clearly 
located in a vicinity of the Stark saddle point. From a 
quadratic approximation of the potential around the sad- 
dle point one expects a harmonic oscillator spectrum of 
energy levels and wave functions. We even find nodal 
patterns in the exact quantum states which correspond 
to "excited" states of the quadratic approximation, and 
thus demonstrate the usefulness of the simple approxi- 
mation and the existence of the resonances predicted by 
a classical treatment of the system 0, 0] • 

In Sec.|TT]we introduce the system as well as the second- 
and third-order power-series expansions of the potential 
in the vicinity of the Stark saddle point. Position space 
representations of exact quantum wave functions which 
prove the existence of resonances located closely at the 
position of the saddle are presented in Sec. Mil Conclu- 
sions are drawn in Sec. IIV1 
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II. HAMILTONIAN AND APPROXIMATIONS 
IN THE VICINITY OF THE STARK SADDLE 
POINT 

A. Hamiltonian and exact quantum calculations 

In Hartree units the Hamiltonian of a hydrogen atom in 
crossed external fields with an electric field / orientated 
along the x-axis and a magnetic field along the z-axis 
represented by the vector potential A has the form 



H=\(p + Af-- + fx. 
2 r 



(1) 



The parity with respect to reflections at the (z = 0)-plane 
is a good quantum number and allows considering states 
with even and odd z-parity separately. 

It is a common and efficient method to rewrite the 
Schrodinger equation for numerical calculations in di- 
lated semiparabolic coordinates [HI, [l3[ , 
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and the generalized eigenvalues A = — (1 + 2b 4 E), which 
are related to the energies E of the quantum states. The 
calculation of the resonances is done with a diagonaliza- 
tion of a matrix representation of the Schrodinger equa- 
tion © and the complex rotation method [14Hlq |. The 
necessary complex scaling of the coordinates r is intro- 
duced via the complex convergence parameter 

b=\b\e li} / 2 . (5) 

Resonances appear as complex eigenvalues E, where the 
real part of E represents the energy and the imaginary 
part is related to the width T = — 2lm(E). 

An adequate complete basis for the matrix representa- 
tion of the Schrodinger equation (PJ is given by 

\n^,n v ,m) = \n^,m) ®\n v ,m), (6) 

where \n e m) are the eigenstates of the two-dimensional 
harmonic oscillator. The position space representation in 
dilated semiparabolic coordinates has the form 



m\)/2]\[{n u -\m\)/2}\ 




m\)/2]\ [{n v + \m\)/2]\ 
(p)f nvm (vy m * (7a) 



with the "radial" wave functions 

f nm {Q)=e-^/ 2 e^L^_ lml)/2 ( Q 2 ) (7b) 

expressed in terms of Laguerre polynomials L"(x). Due 
to the complex scaling /i and v become complex coor- 
dinates. In properly complex conjugated wave functions 
obtained with the complex rotation method one has to 
bear in mind that only the intrinsically complex parts 
have to be conjugated (THE?]]; i.e., in the complex con- 
jugated version ipn ll n v m 01 the wave function (|7af e lmv 
is replaced with e~ lmv and the "radial" parts f n „m{p), 
fn v mky)i which become complex only by the complex 
scaling, are not conjugated. 

It is very effective to perform the matrix setup of the 
Schrodinger equation §5§ completely algebraically by ex- 
pressing Eq. ([3]) in terms of harmonic oscillator creation 
and annihilation operators. Consequently, the expansion 
coefficients Ci n ^ n ^ m of the eigenstates obtained in a ma- 
trix diagonalization 

(fj,,i>, <p), (8a) 

rifj, ,n u ,m 

®i(V,V, ( P)= Cin^mlPn^numi^ViV) ( 8b ) 

,n„ ,m 

belong to position space wave functions (|7al) . which are 
normalized eigenstates of two coupled two-dimensional 
harmonic oscillators, 

oo oo In 
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— &n tl n'i 1 &n v n' 1J &rnm' l (9) 

but are not orthogonal and not normalized in the physical 
position space r = (x, y, z) T . The eigenstates (|8a|) and 
(|8b|) of the Schrodinger equation ([3]), however, can be 
normalized in such a way that 



d d r^(/x,i/, (p)^j{ji,u t ip) 

oo oo 2tt 

= b 6 J dp J dvj dip nv{n 2 + v 2 ) = S i:j . (10) 





Since the proper complex conjugation of the wave func- 
tions effects only intrinsically complex parts, their square 
moduli must be replaced with as introduced in Eqs. 

(f5aj) and (|8b[) to obtain a measure for the probability 
density. For wave functions of decaying states this prod- 
uct will be complex, and we will visualize the modulus 
\tPn ll n v m' l Pn fi n„m\ in Sec. IHII to show where the resonance 
wave functions are located. 
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(a) V(5,y,0) [a.u.] 
0.05 r 




FIG. 1: (Color online) Potential of the hydrogen atom in 
an external electric field (dashed red grids in (a) and (b)) 
in comparison with the second-order (a) and third-order (b) 
power-series expansions around the saddle in the (£,y)-plane 
with £ = x — x s . Both approximations reproduce the potential 
the vicinity of the saddle at the origin correctly, however, both 
are only valid very close to the saddle point. 

B. Approximations in the vicinity of the saddle 
point 

To investigate resonances located at the saddle point, 
the electric potential of the combined contributions of the 
nucleus and the external electric field, 

V f = ~+fx, (11) 

is expanded into a series up to third-order. As was 
shown in Ref. [Hj], already the second-order approxima- 
tion leads to a good one-to-one correspondence with ener- 
gies of exact resonances of the full Hamiltonian (fTJ). Here, 
we additionally consider an approximation including all 
third-order terms to estimate the quality of the approxi- 
mation carried out in the vicinity of the saddle point by 
comparing the second- and third-order results. The sad- 
dle point has the coordinates r s = ( — 1/v 7 /, 0, 0) T and 
its energy has the value Vf(r s ) — —2y/f. Using the co- 
ordinate shift £ = x — x s a power-series expansion around 
r s yields 

-/ 2 e 3 + ^/W + - 2 ) + 0((r-r s ) 4 ). (12) 

In Fig.[T]the second- and third-order power-series expan- 
sions around the saddle point and the full electric po- 
tential (|lip are compared in the (£,y)-plane. The saddle 



structure is clearly visible and the figure demonstrates 
that both approximations are only valid very close to the 
saddle point, and thus can describe resonances correctly 
only in its vicinity. 

The approximated Hamiltonian close to r s reads 

H = - (p\ + p 2 y + pl) - jyp 6 + i 7 2 y 2 

-/ 2 £ 3 + ^/W + ^ 2 ) + 0((r-r s ) 4 ), (13) 

where the gauge A = (—7?/, 0, 0) was used as in Refs. 
05 0> 03 ' since it leads to a simple structure of the terms. 
Introducing complex scaling parameters Sj via 

£ = Qf/sfi, PC = (14a) 
V = Q v /sy, p y = s y P y , (14b) 
z = Q z /s z , p z = s z P z (14c) 

we calculate the resonances of the Hamiltonian with the 
complex rotation method. To do so, a matrix represen- 
tation with a basis of products 

\n£,n y ,n z ) = \n^) ® \n v ) ® \n z ) (15) 

of one-dimensional harmonic oscillator states |rij) is built 
up and diagonalized. 

In the much simpler quadratic approximation, i.e., 
in the case in which the two third-order terms V3 = 
-/£ 3 + 3/ 2 £(y 2 + z 2 )/2 in Eq. Q2) are neglected, the 
approximated Hamiltonian yields a harmonic oscillator 
spectrum 




where all frequencies Ui can be calculated analytically 
[US EH- ^ includes one inverted oscillator, which leads 
to a purely imaginary frequency u)%, and thus represents 
the resonance character of the states. The frequency u) z 
belongs to the z motion of the electron, whereas u\ and 
LO2 belong to new variables introduced via a canonical 
transformation [7j to separate the x and y motions. 

III. RESULTS AND DISCUSSION 

In Ref. [IH it was shown that the energies of some of 
the resonances calculated in the quadratic approximation 
around the saddle point agree very well with the exact 
quantum energies over a large range in the parameter 
space. The region considered was defined by lines, i.e., 
one-dimensional objects, in the two-dimensional param- 
eter space to obtain clearly identifiable results. In this 
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FIG. 2: (Color online) Comparison of the resonance energies 
(real parts of the complex energy eigenvalues) obtained in ex- 
act solutions of the full Hamiltonian (fl]) (red dotted lines) 
and the second- (solid blue lines) and third-order (filled green 
squares) approximations. The results are shown for even (a) 
and odd (b) z parity separately. Position space represen- 
tations of the exact quantum resonances labeled by capital 
letters are shown in Figs. [3] to [5] In (a) and (b) the dashed 
grey lines represent the saddle point energy. 



paper we chose one of these lines, viz. 

7 = 0.008 x a, 
f = 0.0003 x a, 
< a < 1, 



(17a) 
(17b) 
(17c) 



on which all parameters used in what follows are located. 
Figure [2] shows how the real parts of the resonance ener- 
gies behave on this line. The red dotted lines represent 
the exact quantum solution of the full Hamiltonian (JTJ, 
the solid blue lines and the filled green squares denote 
the second- and third-order energies, respectively, and 
the saddle point energy is marked by the dashed grey 
lines. Results are shown for even [Fig. Ufa)] and odd 
[Fig. Mh)} z parity. 

As was already pointed out in Ref. [12| some of the 
second-order resonances are traced by the exact solu- 
tions of (pj . Here we see that this is in particular true for 
all resonances of the second-order power-series expansion 
(solid blue lines in Fig. [5]) , which have a clear correspon- 
dence in the third-order results (filled green squares in 
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FIG. 3: (Color online) Density plots of the resonance labeled 
Al in Fig. [2] in the planes z = 0, y — 0, and x — x s . The 
solid magenta lines mark the saddle point. One can clearly 
see that the probability density of the resonance is restricted 
to a close vicinity around the saddle point. 



Fig. [2]), i.e., which have almost the same energies in both 
approximations for all parameters a. In other words, ap- 
proximated resonances in the vicinity of the Stark saddle 
point whose energies seem already to be converged in 
the two lowest-order power-series expansions have also 
a counterpart in the exact spectrum of the full quan- 
tum system described by the Hamiltonian (JXJ) . In these 
cases the simple power-series expansions seem to pro- 
vide a good approximation for resonances located at the 
saddle. Due to the very good agreement of the second- 
order and the exact quantum results in a large region 
of the parameter space it was already concluded in Ref. 
[IH that the corresponding exact crossed-fields hydrogen 
atom resonances must be closely related with localized 
states predicted by Clark et al. 

In this paper we show using the position space repre- 
sentation of the exact quantum wave functions that they 
are clearly centered at the Stark saddle point in all cases, 
in which the good agreement in the energies described 
above appears. The best correspondence is expected for 
the lowest transition state resonance with quantum num- 
bers n z = 0, n\ = in the second-order approximation 
and for the largest field strengths. The position space 
wave function of this resonance for a = 1, i.e., 7 = 0.008, 
/ = 0.0003 (labeled Al in Fig. [2j is shown in Fig.[3]in the 
planes z = 0, y = 0, and x — x s . All three sections show 
that the probability density of the resonance is centered 
at and restricted to a close vicinity of the saddle point, 
which is marked by the solid magenta lines. The trans- 
formation of the quantum resonances corresponding to 
the same second-order line for decreasing field strengths 
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FIG. 4: (Color online) Density plots of the resonances labeled 
A2, A3, A4, and A5 in Fig.[2]m the z = plane. For decreas- 
ing field strengths one can observe that the resonances are 
decreasingly centered at the saddle point. 



can be observed in Fig. where examples for a — 0.5, 
a = 0.4, a = 0.3, and a — 0.2 are drawn in the z — 
plane. Between a = 1 (resonance labeled Al in Fig. 
12]) and a — 0.5 (A2) the shape of the resonance does 
not change much, it stays centered at the saddle point, 
however, a radially oscillating probability density at the 
borders becomes more and more pronounced indicating 
an increasing contribution of typical Rydberg states. For 
a = 0.4 (A3) one can already observe that the whole wave 
function is affected by the oscillations. The influence of 
the Coulomb potential becomes large enough to lead to 
a wave function with nodal lines in radial direction for 
a = 0.3 (A4), whose main maximum is still located at 
the Stark saddle point. This localization is lost if one 
decreases the field strengths further, which can be seen 
for a — 0.2 (A5), where the maximum of the probability 
amplitude is far away from the saddle. 

There are a number of "excited" second-order states 
which have a correspondence in the third-order power- 
series expansion as well as in the exact quantum spec- 
trum. Figure [5] compares resonances with excitations in 
the n z quantum number. In the second-order approx- 
imation these energies correspond to excited harmonic 
oscillator states in z direction. The resonance labeled B 
represents the quantum number n z = 1, i.e., the first ex- 
cited state. One nodal line at z = is expected in the sec- 
tion and can be observed in Fig. [5] One may argue that a 
nodal plane at z = is not surprising for a state with odd 
z parity, however, further excitations whose nodal pat- 
terns reproduce the expectations correctly can be found. 
Comparing the second-order power-series result for quan- 
tum numbers n z = 2, ri\ = and its third-order counter- 
part in Fig. (5Ja) we find two possibilities to assign exact 
quantum resonances. Indeed, we find the two states la- 
beled CI and C2 whose probability densities are located 
very close to the saddle point. Both of these resonances 
show the two nodal lines expected for n z = 2 and demon- 



FIG. 5: (Color online) Density plots of the resonances B, CI, 
C2, and D (cf. Fig.[2| in the z = plane. The resonances can 
be assigned to second-order approximations with n\ = and 
n z 7^ 0. The nodal patterns expected for excited states in z 
direction can be observed. 



strate that a rich spectrum of states dominated by the 
local shape of the potential around the Stark saddle point 
is present. The resonances CI and C2 have almost the 
same real part of the energy but differ significantly in 
their imaginary parts. One can assume that these reso- 
nances can be assigned to two second-order states with 
n z = 2, n\ = and different ni values, which do not 
influence the real energy in the second-order approxima- 
tion. A definite answer, however, is not possible because 
the excellent agreement in the real parts of the energies of 
the approximated states with the exact quantum results 
of the full Hamiltonian (|T|) does not hold for the imag- 
inary parts 12 1 . The imaginary parts of the energies 



obtained with the second- and third-order power-series 
expansions of the potential do not agree as well as their 
real parts and cannot be considered converged in these 
simple approximations. Their quality does not allow for 
a comparison with their numerically exact counterparts. 
Note that the exact computations with the complex ro- 
tation method provide accurate results for both the real 
and imaginary parts of the energies. Finally one can even 
observe the resonance labeled D with three nodal lines as 
one would expect for a n z = 3 state, and indeed, the res- 
onance energy can be mapped to that of the second-order 
state n z = 3, n\ = 0. The line belonging to resonance 
n z = 4, n\ = in Fig. [2ja) does not seem to have a 
correspondence in the third-order approximation and it 
was not possible to find exact quantum resonances with 
the proper nodal structure. 

Further examples including excitations rt\ ^ can 
be found in Fig. [5] Due to the energy diagram in Fig. 
[5{a) resonance E can be assigned to quantum numbers 
n z = 0, rii = 1. The position space representation in 
Fig.|5]reveals a complicated shape of the probability den- 
sity, which is not very surprising as the excitation in n\ 
belongs to a coordinate whose origin is in a canonical 
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FIG. 6: (Color online) Density plots of the resonances E, F, 
and G (cf. Fig. [2]) corresponding to states of the quadratic 
approximation excited in m . Complicated nodal patterns are 
found. 



transformation including the positions x, y as well as the 
momenta p x and p y . However, the resonance is clearly 
located close to the position of the saddle. This is also 
true for resonance F, which is assigned to the second- 
order quantum numbers n z = 0, ri\ — 2 and shows even 
a more complicated nodal pattern in position space. Res- 
onance G is connected with an excitation in both n z and 
n\ (n z = 1, n\ = 1) and exhibits a nodal plane at z = 
due to the respective quantum number n z — 1. 



IV. CONCLUSION 

By determining the position space representation of 
quantum mechanically exact wave functions we have 
proved the existence of resonances located in a close 
vicinity of the Stark saddle point which were predicted 
by simple classical approximations as, e.g., the transition 
state theory. The results show that near the saddle one 
can find a spectrum of several states belonging to the 
local neighborhood of the saddle and ignoring the global 
structure of the Coulomb potential. It reproduces the 



expectations of a simple quadratic power-series expan- 
sion around the saddle in the energies as well as in the 
shape of the wave functions. In particular, one can, in 
the quadratic approximation, separate three one dimen- 
sional harmonic oscillators one of which is connected with 
the z direction. In this spatial direction we found nodal 
planes of resonances which have the structure of ground 
and excited states from n z = up to n z = 3. Com- 
plicated nodal patterns are found for excitations in the 
quantum number ri\ not directly connected to a spatial 
coordinate. 

The comparison of the two approximations of the reso- 
nance energy eigenvalues (cf. Fig. [2]) showed that some of 
the second-order energies have a distinct correspondence 
in the third-order approximation while some have not. 
In all cases in which an assignment of the results of both 
approximations is possible we were also able to detect 
an exact quantum resonance with a strong localization 
at the saddle, whereas such a strong connection to the 
saddle could not be observed for exact resonances not be- 
longing to a pair of second- and third-order resonances. 
Despite this clear result it must be noted that, as was 



already pointed out in a previous publication 12j, the 
simple power-series expansions are not capable of repro- 
ducing the resonance widths (or imaginary parts of the 
energy eigenvalues) correctly. We have demonstrated, 
however, that a convergence of the approximated reso- 
nance energies (or real parts) is already a strong signa- 
ture of the existence of exact quantum states centered at 
the saddle. 

The results prove that resonances located at an outer 
saddle predicted in a variety of theoretical work 0, 0, 0, 
H, 0] and supported by experimental results @ do exist 
in the Coulomb potential superimposed by two external 
fields. The correspondence already appears for very sim- 
ple second- and third-order power-series expansions of 
the potential at the saddle point. For future work com- 
parisons with more thorough approximations applicable 
to the problem, e.g., the normal form expansion devel- 
oped for identifying the classical transition state in high- 
dimensional systems [7|, [lfj, or its quantum analogue 
[H E3 will be of high value. 



[1] S. K. Bhattacharya and A. R. P. Rau, Phys. Rev. A 26, 
2315 (1982). 

[2] J. C. Gay, L. R. Pendrill, and B. Cagnac, Phys. Lett. A 

72, 315 (1979). 
[3] A. R. P. Rau and L. Zhang, Phys. Rev. A 42, 6342 (1990). 
[4] C. W. Clark, E. Korevaar, and M. G. Littman, Phys. 

Rev. Lett. 54, 320 (1985). 
[5] M. Fauth, H. Walther, and E. Werner, Z. Phys. D 7, 293 

(1987). 

[6] G. Raithel, M. Fauth, and H. Walther, Phys. Rev. A 47, 
419 (1993). 

[7] T. Uzer, C. Jaffe, J. Palacian, P. Yanguas, and S. Wig- 



gins, Nonlinearity 15, 957 (2002). 
[8] C. Jaffe, D. Farrelly, and T. Uzer, Phys. Rev. Lett. 84, 
610 (2000). 

[9] C. Jaffe, D. Farrelly, and T. Uzer, Phys. Rev. A 60, 3833 
(1999). 

[10] S. Wiggins, L. Wiesenfeld, C. Jaffe, and T. Uzer, Phys. 

Rev. Lett. 86, 5478 (2001). 
[11] H. Waalkens, R. Schubert, and S. Wiggins, Nonlinearity 

21, Rl (2008). 

[12] H. Cartarius, J. Main, and G. Wunner, Phys. Rev. A 79, 

033412 (2009). 
[13] J. Main and G. Wunner, J. Phys. B 27, 2835 (1994). 



7 



[14] W. P. Reinhardt, Ann. Rev. Phys. Chem. 33, 223 (1982). 
[15] D. Delande, A. Bommier, and J. C. Gay, Phys. Rev. Lett. 

66, 141 (1991). 
[16] N. Moiseyev, Phys. Rep. 302, 212 (1998). 
[17] T. N. Rescigno and V. McKoy, Phys. Rev. A 12, 522 



(1975). 

[18] R. Schubert, H. Waalkens, and S. Wiggins, Phys. Rev. 
Lett. 96, 218302 (2006). 



